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A  GENERAL  MATHEMATICAL  MODEL  FOR  BEAM  AND  PLATE  VIBRATION 
IN  BENDING  MODES  USING  LUMPED  PARAMETERS 

SUMMARY 


[The  analytic  method  proposed  in  this  paper  solves  for  natural  frequencies  and 
mode  shapes  of  distributed  mass  and  stiffness  systems  such  as  beams  and  plates.  The 
structure  is  divided  into  an  equivalent  system  of  discrete  masses  and  springs.  A  sche¬ 
matic  diagram  can  be  built  up  to  show  the  nature  of  each  item  affecting  the  system. 
Increased  off-hand  evaluation  of  design  changes  results.  A  group  of  simultaneous 
equations  is  made  up  from  the  discrete  schematic  system.  These  equations  comprise 
the  mathematical  model  of  the  systen^ 

The  mathematical  model  of  the  beam  is  completely  general.  The  model  is  a 
matrix  equation  which  will  accept  any  input,  any  restraint,  or  any  combination  of  end 
conditions.  One  model  will  suffice  for  all  beam  problems.  Damping  can  be  considered, 
but  is  not  included  in  the  main  thesis  of  the  paper.  The  model  is  in  sufficient  detail  to 
accept  the  particulars  of  a  nonuntform  beam  where  mass  and  stiffness  can  be  broken 
into  lumped  discrete  but  unequal  values.  The  number  of  masses  limits  the  eigenvalue 
solutions.  Reasonable  accuracy  is  dependent  on  a  minimum  of  three  masses  between  nodes. 
This  is  equivalent  to  four  masses  per  modal  number,  plus  one  for  the  end  node  .  Thus 
a  system  with  seventeen  masses  would  be  reasonable  for  the  calculation  of  the  first  four 
modes  of  a  beam. 

Flat  plates  are  simulated  as  two  sets  of  cross-coupled  beams.  A  nonuniform 
plate  model  is  developed.  End  conditions  of  the  individual  beams  may  be  varied,  masses 
may  be  equated  to  zero,  etc,  to  set  up  a  rather  complicated  model.  As  in  the  beam  model, 
the  plate  model  is  very  general  and  only  one  model  is  necessary  where  the  size  is  large 
enough  to  accept  the  required  detail.  Matrix  size  limits  the  number  of  eigenvalue 
solutions  available.  The  model  of  a  cross-coupled  plate  of  sixteen  masses  is  derived  in 
the  paper.  Actually,  practical  solutions  would  require  a  much  larger  equation. 

{jhe  model  does  not  account  for  any  extensional  deformation.  Only  bending  is 
considered.  The  application  of  the  method  to  shells  and  other  curved  plates  is  suggested 
when  only  bending  modes  need  to  be  considered 


SECTION  I.  INTRODUCTION 


Mechanical  systems  comprised  of  distributed  mass  and  distributed  stiffness  have 
an  infinite  number  of  normal  modes  of  vibration.  Beams  are  characteristic  of  the  distri¬ 
buted  systems  and  have  only  one  set  of  normal  modes  per  single  direction  of  deflection. 
They  therefore  are  representative  of  the  simplest  approach  to  distributed  system  dynamic 
analysis.  Other  systems  such  as  plates  and  shells  can  be  represented  by  a  system  of 
cross-coupled  beams  once  a  suitable  system  of  impedance  and  coordinate  designation  is 
selected. 

Systems  comprised  of  uniformly  distributed  mass  and  stiffness  lend  themselves 
to  the  usual  methods  of  solution,  by  utilizing  some  form  of  the  Rayleigh  method  by  equa¬ 
ting  maximum  strain  to  maximum  kinetic  energy  for  each  mode.  These  methods  arrive  at 
mode  shapes  and  natural  frequencies  but  not  finite  response  amplitudes.  The  incorpora¬ 
tion  of  damping  into  the  calculations,  then,  allows  an  approach  to  finite  amplitudes  of 
response  under  steady-state  conditions  by  equating  input  energy  to  dissipated  energy. 

Some  systems  are  not  truly  uniformly  distributed  but  can  be  simulated  by  an  equivalent 
uniformly  distributed  system.  Typical  of  this  class  is  a  plate  with  equally  spaced  like 
stiffeners.  The  lower  numbered  modes  of  this  type  of  stiffened  plates  will  respond  in 
close  approximation  to  mode  shapes  of  uniformly  constructed  plates. 

Some  structures  have  irregular  details  of  mass,  stiffness,  or  end  conditions  that 
do  not  allow  a  convenient  method  of  solution  by  distributed  system  analysis .  They  must 
be  simulated  by  a  system  of  lumped  masses  influenced  by  local  spring  effects.  The  use 
of  these  lumped  methods  is  usually  limited  to  a  few  of  the  lowest  natural  frequencies 
because  a  lumped  system  has  a  finite  number  of  degrees  of  freedom  and  the  deflection  of 
each  mass  reflects  a  single  point  on  the  mode  shape  curve.  Good  definition  of  mode  shape 
requires  a  minimum  of  about  four  masses  between  modes.  Practical  considerations  of 
the  problem  solutions  limits  the  number  of  useable  lumped  masses.  Since  design  prob¬ 
lems  are  usually  associated  with  low  numbered  modal  response,  the  lumped-mass 
method  could  be  of  considerable  importance. 

Any  solution  of  a  distributed  system  dynamic  problem  is  of  a  complicated  nature. 

A  method  that  allows  a  simple  subjective  comprehension  of  each  of  the  details  has  con¬ 
siderable  value  to  the  analyst.  It  will  increase  off-hand  evaluation  of  the  effect  of  any 
change,  and  lead  to  more  effective  design  considerations.  The  analytic  method  proposed 
in  this  paper  takes  each  specific  detail  and  transforms  it  into  an  easily  comprehendable 
part  of  the  system  diagram.  The  method  is  set  up  so  that  it  is  completely  general.  One 
program  could  be  developed  to  accept  all  beam  problems.  Similarly,  one  program  could 
handle  all  plate  solutions. 

The  method  will  be  derived  by  detailed  steps  so  that  this  paper  can  serve  as  a 
ready  reference  to  users  of  dynamic  analysis. 
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SECTION  n.  METHOD  DERIVATION 


The  first  step  in  the  method  derivation  is  to  review  the  basic  equations  relating 
moment  to  deflection  (Fig.  1).  When  the  deflection  is  very  small  and  there  is  no  net 


FIGURE  1.  SECTION  OF  A  BENT  BEAM 

elongation  of  the  beam,  the  deflection  of  the  surface  of  the  beam  (BB',  Fig.  1)  can  be 
considered  parallel  to  the  undeflected  axis  of  the  beam,  and  the  deflection  of  the  neutral 
axis  of  the  beam  (00')  may  be  considered  perpendicular  to  the  axis  of  the  beam.  If  we 
arbitrarily  set  AL  equal  to  OB,  then  simple  geometric  analysis  will  show  that  BB'  equals 
00'  in  amplitude,  and  that  they  are  90  degrees  apart  in  direction.  The  equation  for  the 
moment  acting  on  section  B'B'  is 


■M  = 


El 


Ae 

AL 


(1) 


where 


moment 

E  =  Young's  Modulus 
BB* 

A©=  — T"  in  radians  for  this  case 
AL 

I  =  moment  of  inertia 
AL  =  unit  length  of  beam. 

Now  the  moment  is  equal  to  the  product  of  force,  F,  times  arm  length,  AL.  Equation 
( 1 )  can  be  written 
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(2) 


F  AL  =  E  I 


BB' 

(AL)2 


E  I 


00' 

(AL)2 
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or 


F  = 


00' 


(AL) 


^E  I. 


(2a) 


Since  the  stiffness  of  the  beam,  K,  is  equal  to  F  divided  by  the  deflection  00', 
F  El 


K  = 


00 


(AL) 


(3) 


Equation  (3)  has  the  restrictions  that  the  deflection  be  small  and  that  moments  acting  on 
the  stiffness  must  be  acting  through  the  arm  length  of  AL.  Figure  I  can  now  be  simula¬ 
ted  by  a  discriptive  model  (Fig.  2) . 


FIGURE  2.  BEAM  ELEMENT  MODEL 


To  simulate  a  beam.  Figure  2  can  be  modified  to  represent  a  series  of  linked  massless 
trusses  with  stiffness  and  mass  lumped  arbitrarily  and  with  the  basic  concepts  of  Figure 
2  adhered  to.  Such  a  system  is  shown  in  Figure  3.  It  must  be  remembered  that  all 
effective  motion  of  the  spring  attachment  points  is  in  the  horizontal  direction,  and  that 
of  the  hinges  in  the  vertical  direction  (each  related  to  the  other) .  The  interrelationship 


K  =  Spring  constant 
M  =  Mass 
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FIGURE  3.  SIMULATED  BEAM 


of  these  two  displacements  leads  to  a  defining  of  the  spring  deflection  in  terms  of  the 
deflection  in  the  vertical  axis.  Each  truss  is  hinged  at  the  base  to  the  adjacent  trusses. 
As  defined  above,  this  hinge  allows  motion  only  in  the  vertical  direction.  Each  hinge 
deflection  is  representative  of  the  displacement  of  the  neutral  axis  of  the  beam  at  its 
relative  location.  The  relative  displacement  of  one  hinge  of  a  truss  to  the  other  is  equal 
in  magnitude  to  the  horizontal  displacement  of  the  spring  attachment  point.  Thus  all 
spring  deflections  can  be  equated  to  actual  beam  deflections.  Let  us  define  the  net 


r 

FIGURE  4.  SPRING  DEFLECTION  MODEL 

deflection  of  spring  A  B,  Figure  4,  by  and  the  displacement  of  A  by  A^,  etc. 

The  following  equation  develops 


Aa  =  A2  -  Ai, 

Ag  =  A3  -  Aj, 

A,^  =  A^  -  A  =  A3  -  2A2  +  Ai.  (4) 

AB  B  A  ^  '  ' 


The  spring  constants  and  the  masses  do  not  have  to  be  equal  for  analytic  solution. 
However  the  solution  calculations  are  greatly  simplified  when  all  trusses  are  of  equal 
size,  except  for  the  halves  at  the  ends.  This  in  effect  divides  the  beam  into  segments 
of  equal  length  between  hinges,  with  the  ends  being  half  segments.  Figure  5  depicts  the 


FIGURE  5. 


(B)  HINGED  END 


(C)  FIXED  END 
END  CONDITIONS 
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three  types  of  end  conditions  usually  assumed  in  beam  problems.  The  graphic  descrip¬ 
tion  is  self  evident.  Now  comes  some  of  the  real  value  of  this  method.  Most  deviations 
from  the  standard  end  conditions  can  be  simply  worked  into  the  solution.  Figure  6  shows 
a  few  of  the  possible  non-standard  end  conditions  of  practical  consideration.  The  non¬ 
standard  impedance  values  attached  to  the  ends  are  easily  added  to  the  general  equations 
so  that  this  type  of  problem  is  included  in  the  general  solution. 


(D)  RESTRAINED  HINGE 


(E)  FLEXIBLE  ATTACHMENT 


(F)  END  MASS  RESTRAINT 


FIGURE  6.  NON-STANDARD  END  CONDITIONS 


Attachments  to  the  beam  also  can  be  accounted  for  in  the  general  solution.  Fig¬ 
ure  7  shows  the  diagramatic  inclusion  of  such  attachments  to  the  simulated  system. 


(G)  FIXED  ATTACHMENT  (H)  SPRING  ATTACHMENT  (I)  MASS  ATTACHMENT 


FIGURE  7.  BEAM  ATTACHMENTS 


It  should  be  noted  that  damping  has  not  yet  been  inserted  into  the  diagram.  A 
type  of  damping  simulation  is  possible  but  not  very  practical  within  the  present  state-of- 
the-art.  Its  usefulness  will  depend  upon  further  research.  The  diagramatic  model  de¬ 
scribed  up  to  this  point  is  for  free  vibration  systems.  The  characteristics  of  the  system 
can  be  conveniently  represented  by  impedance  values.  Such  representation  can  lead  to 
an  understanding  of  the  damping  properties  of  the  system.  Structural  damping  is  inserted 
into  the  impedance  equation  as  follows; 

Z  =  ju)M 


=  jwM  +^. 

JO)  m 


(5) 


M  =  mass 
K  =  stiffness 
7]  =  damping  loss  factor 
ct»  =  angular  frequency. 

Force,  F,  is  related  to  displacement,  A,  by  the  relation 

F  =  (Z)  (ja;A).  (6) 

The  simulated  system  shown  in  Figure  3  can  be  solved  by  using  the  usual  equations  of  con¬ 
tinuity  of  force  and  moment.  Damping  can  be  included  in  the  impedance  equation.  How¬ 
ever  when  used,  it  is  recommended  to  solve  the  reactive  portion  first  without  damping, 
and  then  insert  damping  to  calculate  amplitudes.  For  simple  beam  and  plate  models, 
damping  can  be  schematically  added  as  in  Figure  8.  Any  of  the  usual  damping  concepts 
are  useable  but  here  we  use  a  damping  symbol  that  is  indicative  of  slipping  or  Coulombs^ 


damping  and  this  type  is  proportionate  to  displacement.  Thus  the  damping  is  proportionate 
to  the  same  displacement  that  works  on  the  spring  located  between  each  of  the  trusses. 


The  loss  factor  is  unit  length  or  area  dependent.  Once  the  mode  shape  is  determined 
for  a  mode  the  input  energy  must  equal  the  output  energy, 


(7) 


where  m  is  representative  of  forced  input  locations  and  n  is  representative  of  response 
deflections  of  the  various  springs.  Equation  (7)  would  be  applicable  to  beams  or  plates 
without  attachments.  By  including  representative  attachment  loss  factor  increments, 
equation  (7)  could  become 


s 


jTi  K  a2  + 
■'  '  n 


N' 


(8) 


where  N',  77'  and  A^,  refer  to  the  damping  caused  by  the  attachments.  The  schematic 
inclusion  of  the  damping  caused  by  the  attachments  is  represented  in  Figure  9.  Equation 
(8)  will  solve  for  discrete  values  of  amplitude  for  each  mode.  The  A  deflections  are 
dependent  on  the  y  displacements.  This  aspect  will  be  discussed  after  development  of 
the  general  solution  equation. 


FIGURE  9.  INCLUSION  OF  ATTACHMENT  DAMPING 

It  has  been  the  purpose  here  to  show  the  possibility  of  damping  analysis  by  using 
the  analytical  model.  The  remainder  of  the  report  will  deal  with  undamped  systems  only. 


SECTION  m.  BEAM  SOLUTIONS 


Development  of  a  General  Beam  Matrix  Equation 

The  beam  solution  uses  a  set  of  equations  that  define  the  continuity  of  force  and  of 
moments.  Figure  10  is  a  typical  model  of  a  beam  of  six  elements  of  mass.  All  motion 
of  the  beam  is  assumed  in  the  ±y  directions  only.  End  restraints  may  impose  a  reaction 
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FREE 

FREE 

^0  = 

=  0,  Ro  -  0 

Ivly  ”0,  Ry  —  0 

FIXED 

FIXED 

yo  = 

0 

II 

0 

ye  =  0.  y?  =  0 

HINGED 

HINGED 

yo  = 

0,  -M])  =  0 

yy  =  0,  -My  =  0 

Mq,  ^7  are  end  moments,  in.  lbs. 

Rfl,  R7  are  end  support  forces,  lbs. 

Kj  etc.  is  stiffness  constant  =  jf ,  lbs. /in. 

Aj  etc.  is  the  net  deflection  of  the  spring  "K",  in. 
yo  etc.  is  the  neutral  axis  deflection  of  the  beam,  in. 
i  is  the  height  and  base  of  each  full  triangle ,  in. 

Mj  etc,  is  the  lumped  mass  of  the  section  of  the  beam. 


force  R  and/or  a  moment  of-M:  The  end  may  have  a  displacement  of  yg  or  yy.  Each 
element  of  mass  M  and  stiffness  K  may  have  its  own  assigned  quantitative  value.  In  the 
case  of  a  uniformly  distributed  system,  all  M  values  would  be  equal  and  all  K  values 
would  be  equal. 


The  Symbols  used  in  the  derivation  of  the  equations  are  defined  below  Figure  10. 
The  counter-clockwise  moment -ftto  is  considered  positive,  whereas  the  clockwise  mo¬ 
ment  My  is  considered  positive.  This  is  the  usual  convention  used  in  beam  equations. 
The  deflection  of  y  in  the  downward  direction  is  considered  positive  and  the  deflections 
of  the  springs  to  the  right  are  considered  positive.  Both  support  forces,  Ro  and  Ry,  are 
considered  positive  in  the  upward  direction. 
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The  first  set  of  equations  to  be  evolved  are  based  on  the  reasoning  leading  to  equa¬ 
tion  (4) .  We  define  Aj  as  the  deflection  imposed  on  and  A2  as  the  deflection  imposed 
upon  K2,  etc. ,  a  positive  value  of  A  reflects  an  elongation  of  the  spring.  Remember 


that  the  moment  arm  of  the  springs 

is  equal  in  length  to  the  distance  between  hinges. 

0  =  -Ai  +  y2  -  3yi  +  2yo 

(9) 

0  =  -A2  +  ys  -  2y2  +  yi 

(10) 

0  =  -A3  +  y4  -  2y3  +  y2 

(11) 

0  =  -A4  +  yg  -  2y4  +  y3 

(12) 

0  =  -Ag  +  ye  -  2y5  +  y4 

(13) 

0  =  -Ag  +  2y,  -  Syg  +  yg 

(14) 

In  equations  9  through  14  the  y's  are  actual  deflections  at  the  locations  as  indicated  by 
the  subscript. 


Next, in  Figure  11,  let  us  look  at  the  moments  relative  to  a  section  1,  1  of  the 
beam  passing  through  location  1.  The  deflection  in  spring  Kj  is  induced  by  equal  and 
opposite  moments  on  either  side  of  the  section  1, 1.  For  convenience  let  us  take  the 
moments  to  the  left.  Then  the  continuity  of  moments  at  section  1, 1  will  be: 

Aj  Kj  I  =  M()  -  2  R(). 

This  can  be  resolved  to  a  more  standard  form  of: 

“^0 

0  =  Ai  Ki  -  —  +  .5  Rq.  (15) 

Note  that^o  result  in  moments  of  opposite  sign  as  they  resist  each  other. 

Section  2,2,  Figure  12,  now  involves  more  than  end  conditions  in  the  moment 
equation.  Since  the  springs  at  locations  other  than  the  designated  section  location  only 
transmit  moment,  they  are  not  involved  in  the  equation.  -Mq  and  Rg,  as  external  effects, 
do  not  represent  physical  characteristics  of  the  beam.  The  mass,  on  the  other  hand, 
is  an  impedance  value  and  can  be  treated  as  such.  The  moment  about  section  2,2  in¬ 
duced  by  any  impedance,  Z^,  with  a  net  displacement  of  yj  would  be: 
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FIGURE  11.  BEAM  SECTION  AT  LOCATION  1 


FIGURE  12.  BEAM  SECTION  AT  LOCATION  2 


11 


^(Zl)  =  (Zi)  (j  wyj)  i. 


For  Zi  in  this  case  we  can  substitute  j  co  Mj,  then:  yj  i.  The  moments 

at  section  2,2  can  now  be  equated  as: 

^0 

0  =  A2  Kg  -  —  +  1.  5  Ro  -  0)2  y^  Ml.  ( 16) 

Following  the  same  reasoning  as  above,  equations  of  moment  continuity  can  be  written 
for  sections  at  3,  4,  5,  etc.  For  a  six  mass  beam  the  set  of  moment  equations  would  be: 


0  =  Ai  Ki  -  Y-  + 


5  R, 


(15) 


0  =  A2  K2  -  "T—  +  f.  5  Ro  -  0)2  yi  Ml 


0  -  A3  K3  -  -J—  +  2.  5  Ro  -  2ct)2  yi  Ml  -  0)2  y2  Mg 


0  —  A4  “•  ^  +3.5RQ  —  3CO  Yj  IVI-j  —  2cj  M2  "*  co  yg  Mg 


-M- 

0  =  Kg  -  +  5.  5  Rfl  -  50)2  _  4^2  _  3^2  y^  ^  _  £^2  y^  ^ 

-  0)^5  M5 

We  can  also  write  the  moment  equation  at  y7  and  it  becomes: 


(16) 

(17) 

(18) 


0  =  Ag  K5  -  —  +  4.  5  Rq  -  40)2  yi  Ml  -  3o)2  y2  Mg  -  2o)2  yj  M3  -  0)2  yi  M4  (19) 


(20) 


"^^7  "^0 

0  =  -  ^  +  6  Rfl  -  5.  5  0)2  yi  Ml  -  4.  5  o)2  yg  Mg  -  3.  5  o)2  yj  M3 

(21) 

-  2.  5  0)2  yi  M4  -  1.  5  0)2  yg  Mg  -  0.  5  0)2  yg  Mg, 


If  we  substitute  equation  9  through  14  into  equations  15  through  21  as  applicable, 
we  will  end  with  the  following  list  of  unknowns: 

yo>  drtg,  Rg,  yi,  yg,  yg,  y4,  yg,  yg,  y^j-My,  R^. 
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Any  choice  of  end  conditions  will  reduce  this  group  by  four,  which  gives  eight  unknowns. 
Thus  far  we  have  developed  seven  equations,  numbered  15  through  21.  We  need  one 
more  to  satisfy  our  simultaneous  equation  solution.  In  forced  vibration  the  sum  of  the 
forces  acting  in  the  vertical  direction  must  equal  zero.  The  vertical  forces  are  induced 
by  the  end  forces  Ro  and  R7,  and  the  inertial  forces  by  the  masses  or  any  other  restraint 
at  the  hinges.  The  eighth  equation  can  now  be  developed: 

0  =  -  Ro  -  R7  +  yj  Ml  +  0)2  y2  M2  +  0)2  yj  Mg  +  0)2  y4  M4  +  0)2  yj  M5  +  0)2  yg  Mg 

(22) 


Combining  equations  15  through  22  produces  matrix  equation  (23). 


.5 

0 

0 

+2^1 

-3Ki 

^K, 

0 

0 

0 

0 

1.5 

0 

0 

0 

+  (K2-w2Mi) 

-2K2 

+K2 

0 

0 

0 

0 

2.5 

0 

0 

0 

-2a;2Mi 

+  {K^-C0^M2) 

-2K3 

+K3 

0 

0 

0 

3.5 

0 

0 

0 

-2io^M2 

+(K4-ai2M3) 

-2K4 

+K4 

0 

0 

4.  5 

0 

-1 

0 

0 

-4a)2Mi 

+  (K5-a;2M4) 

-2K5 

+K5 

0 

5.  5 

0 

-f 

0 

0 

-5w2Mi 

-4^2  M4 

-2cj2^4 

+  (K6-Co2M5) 

-3Ke 

+2K6 

6 

0 

+i“ 

0 

-5. 

->4.5co2M4 

-3.  5a;^M3 

-2.5a;2M4 

0 

-1 

-1 

0 

0 

0 

+Cd^M2 

+W^M3 

+C(;^M4 

+cu^Mg 

0 

Matrix  Equation,  Beam  General  Model  of  Six  Masses  Equation  (23) 


Rq 

R? 

Mo 

M7 

yo 

yi 

yz 

ys 

y4 

ys 

ye 

y? 


Equation  (23)  represents  a  completely  general  matnemaucal  modei  of  a  beam. 
The  rectangular  matrix  and  the  right-hand  column  matrix  will  reduce  to  solvable  square 
and  column  matrices  of  the  same  order  upon  selection  of  the  end  conditions.  Note,  from 
the  table  below  Figure  10,  that  each  end  condition  has  two  of  the  possible  variables  equal 
to  zero.  The  selection  of  end  conditions  would  then  reduce  the  rectangular  matrix  of 
equation  (23)  to  a  square  8x8  matrix  and  the  right-hand  column  matrix  will  have  four 
of  its  variables  equal  to  zero  and  therefore  cancel  out  to  the  8th  order. 
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Since  any  selection  of  end  conditions  will  reduce  the  variables  by  two,  the  end 
conditions  need  not  be  alike.  The  matrix  equation  will  accept  any  combination  of  end 
conditions.  This  is  a  distinct  advantage  of  this  type  of  approach.  Once  the  computer  is 
programed  to  solve  the  Eigenvalue  and  Eigenvector  problems  of  full,  non-symmetrical, 
matrices  the  analytical  model  will  produce  solutions  not  convenient  by  other  methods. 
Since  the  matrix  is  not  dependent  on  uniform  constants,  a  non-symmetrical  physical  con¬ 
figuration  may  be  solved.  In  a  formal  sense  of  application  this  method  would  apply  only 
to  stepped  structures  with  uniform  constants  between  steps.  In  a  loose  sense  the  triang¬ 
ular  segments  may  represent  a  lumping  of  the  variable  constants  into  equivalent  values. 
For  instance  if  the  mass  of  the  beam  is  equatable  to  a  function  of  position  on  the  beam. 


dx 


where  =  Mass  of  the  nth  segment  of  the  beam 

p  =  Mass  density  of  the  beam  in  mass  per  unit  length 
X  =  distance  along  the  neutral  axis  of  the  beam. 

The  deflection  of  the  beam  must  be  considered  linear  and  estimates  made  must  be 
for  the  spring  constant  of  each  segment  of  an  irregular  beam. 

Equation  (24)shows  the  reduction  of  equation  (23)  to  represent  a  uniform  symme¬ 
trical  beam  of  six  masses.  This  specific  model  is  used  to  illustrate  the  3  usual  beam 
conditions  of 


1.  Fixed-Fixed 

2.  Hinged-Hinged 

3.  Free-Free 


This  simplified  model,  equation  (24) ,  was  used  because  it 


.5 

-i 

4-2K 

-3K 

+K 

0 

1.5 

-1 

-i 

0 

(K  -  0)%) 

-2K 

+K 

0  = 

2.  5 

0 

-2w% 

+(K  -  w^M) 

-K 

3.  5 

0 

-3c4;2m 

+(K  -  2w%) 

-(K  +  co^M) 

Rfl.T 

^,7 

yo,7 

yi,6 


Matrix  Equation,  Symmetrical  Uniform  Beam  General  Model  of  6 
Masses  Equation  (24) 


y2,5 

y3,4 
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allowed  reasonable  (4  by  4)  hand  solution  to  show  examples  of  problem  solution.  Higher 
order  solutions  are  recommended  for  the  practical  use  of  the  mathematical  model. 

The  solution  of  matrix  equations  such  as  equation  (23)  is  somewhat  involved. 
When  the  rectangular  matrix  is  reduced  to  a  square  matrix  it  represents  the  impedance 
of  the  system.  At  each  resonant  frequency  the  reactive  portion  of  the  impedance  is  equal 
to  zero.  Since  only  reactive  impedance  has  been  inserted  into  the  rectangular  matrix, 
the  values  of  co  that  equate  it  to  zero  are  solutions  of  the  natural  frequencies.  These 
values  of  oo  are  known  as  Eigenvalues. 

The  variables  in  the  column  matrix  represent  system  response.  Each  natural 
frequency  will  have  a  finite  set  of  response  values.  Solutions  of  the  column  matrix  are 
known  as  Eigenvectors.  Even  after  the  insertion  of  a  natural  frequency  into  the  square 
matrix,  equation  (23)  is  not  solvable.  However,  if  we  normalize  the  response  values  to 
one  of  the  Eigenvectors,  e.g.  ,  y^,  then  the  insertion  of  unity  into  the  column 
matrix  makes  the  equation  solvable.  This  means  that  only  relative  displacements  can 
be  determined  from  an  Eigenvalue  -  Eigenvector  solution.  These  relative  values  de¬ 
scribe  a  mode  shape  for  each  natural  frequency. 

The  physical  reason  finite  values  of  response  amplitude  cannot  be  had  from  a 
solution  of  equation  (25)  is  that  response  amplitudes  of  a  steady  state  forced  system 
are  damping  dependent.  Equation  (25)  contains  no  damping  values. 

Equation  (8)  equates  input  energy  to  dissipated  energy.  Each  deflection  y  is 
equated  as 


y  y  , 

n-^y' 


(25) 


Where  y  is  the  normalization  factor.  All  must  be  determined.  Then  the  substitu¬ 
tion  of  equation  (25)  into  equation  (9)  through  (14)  and  the  resultants  into  equation  (8)  will 
lead  to  solutions  of  finite  deflection  values  for  each  mode.  Usually,  at  the  first  few  of  the 
lowest  natural  frequencies  such  as  those  suitable  for  this  method  of  analysis,  these  de¬ 
flections  will  be  representative  of  total  deflections  at  the  natural  frequencies.  This  is 
true  because  modes  other  than  the  resonant  mode  will  have  very  little  amplitude  when 
only  low  frequencies  are  being  considered. 


Examples  -4x4  Matrix  Eigenvalue  Solutions  of  Uniform  Beams 

The  general  solution,  equation  (23),  can  very  often  be  reduced  in  scope.  This 
is  especially  true  when  the  beam  is  symmetrical  and  uniform.  In  fact,  equation  (24)  is 
the  first  step  in  reduction  to  account  for  uniformity.  When  the  beam  is  symmetrical  it 
has  an  even  number  of  masses  and  the  odd  modes  only  are  considered.  Each  deflection 
on  the  right  of  center  should  have  an  equal  deflection  on  the  left  side.  Also  the  moment 
and  reaction  force  at  the  ends  are  equal.  By  substitution  of  one  set  into  the  other, 
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equation  (24)  can  be  reduced  to  half  size.  Equations  26,  27,  and  28  represent  uniform 
symmetrical  beams  of  six  masses  that  reduce  to  a  4  by  4  system.  These  4  by  4  matrices 
are  easily  calculated  by  hand  and  will  be  used  to  check  out  the  accuracy  of  this  analytic 
method  by  general  comparison  with  standard  solutions. 


Matrix  Equations,  Symmetrical  Uniform  Beam,  Specific  CaSes 


1.  Fixed  -  Fixed: 


.5  +K  0 

^■0, 7 

1. 5  -i  -2K  +K 

Mq  7 
> 

2.5  +(K  - 

yz  5 

9 

3.  5  -i  +(K  -  2u;2m)  -(K  +  co^M) 

y3,4 

2.  Hinged  -  Hinged: 


.5 

-3K 

+K 

0 

Rflj 

1.5 

+(K  -  0)%) 

-2K 

+K 

yi,6 

2.  5 

-2aj^M 

+(K  -  co^M) 

-K 

y2,5 

3.5 

-3a)2M 

+(K  -  2a;2M) 

-(K  +  co^M) 

y3.4 

3.  Free  -  Free: 


2K 

-3K 

+K 

0 

yo,7 

0 

(K  -  0)%) 

-2K 

+K 

yi,6 

0 

-2w^M 

+(K  -  w^M) 

-K 

y2,5 

0 

-Sco^M 

+(K  -  2c<;2m) 

-(K  +  w^M) 

y3,4 

(26) 


(27) 


(28) 
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Equation  (29)  is  a  longhand  solution  of  a  four  by  four  matrix.  The  calculation 
will  be  based  upon  this  solution  in  tabular  form.  It  should  be  noted  that  when  A,  B,  C 
or  D  is  equal  to  zero,  a  portion  of  the  tabulation  need  not  be  performed.  For  example, 
D  equals  zero  in  the  solution  of  the  fixed-fixed  beam,  equation  (25) . 


4  by  4  Matrix  Expansion  Model 

A  B  C  D  I 


E  F  G  H 

F  G  H 

E  G  H 

E  F  H 

E  F  G 

=  A 

J  K  L 

-  B 

I  K  L 

+  c 

I  J  L 

-  D 

I  J  K 

I  J  K  L 

N  P  Q 

MP  Q 

MN  Q 

MN  P 

M  N  P  Ql 

=  AFKQ  +  AGLN  +  AHJP  -  AFPL  -  AJGQ  -  ANKH 

-  BEKQ  -  BGLM  -  BHIP  +  BEPL  +  BIGQ  +  BMKH 
+  CEJQ  +  CFLM  +  CHIN  -  CENL  -  CIFQ  -  CMJH 

-  DE JP  -  DFKM  -  DGIN  +  DENK  +  DIFP  +  DMJG 


TABLE  I.  MATRIX  EXPANSION  OF  EQUATION  (26) 


K^ 

Kw^M  co^M^ 

AFKQ  = 

AGLN  = 

.5  (-i“^)(K  -  (x)%)  (-K  -  w%)  =  +  .5  £"^  K^ 

.  5  ( -2K)  ( -K)  (-i  "^)  =  -S. K^ 

-.  5i  “  Vm^ 

AHJP  = 

.  5  (K)  (-jf”^  (K  -  2co2m)  =-.5^"^k2 

+I  ~^w^MK 

-AFPL  = 

-.  5(-i”^)  (K  -  2co2m)  (-K)  =-.5jf“*K2 

+i  co^MK 

-AJGQ  = 

-.  5(-i"^)(-2K)(-K -a)%)  =+i"^K2 

+i"^  co^MK 

-ANKH  = 

-.  5(-jf  “^)  (K  -  0;%)  (K)  =  +  .  K^ 

-  .  5 

-BEKQ= 

-BGLM  = 

+lt~^  (1.5)  (K  -  ci;%)  (-K  -  co^M)  =  -  1.  5i“^  K^ 

+  (-2K)(-K)(3.  5)  =  \-ir^K^ 

+  1.  5i  "^co^M^ 
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Table  I  (Cont'd) 


-BHIP  =  +  i  "Nk)  (2.  5)  (K  -  2aj2]yi) 
+BEPL  =  -i  "^  ( 1.  5)  (K  -  2aj%)( -K) 
+BIGQ  =  -i  2.  5)  ( -2K)  ( -K  -co^M) 
+BMKH  =  -i"^  (3.5)(K  -  w^mXK) 
+CEJQ  =  +K(  1.  5)  ( -i  ( -K  -  co^M) 

+CFLM  =  +K( -a  ( -K)  ( 3.  5) 

+CHIN  =  +K(K)  (2.  5)  (-i 
-CENL  =  -K(  1 .  5)  ( -£  ”‘)  ( -K) 

-CIFQ  = -K(2.5)  (-i"^  (-K  - 

-CMJH=  -K(3.5)  (-i“^  (K) 


=  +  2. 

-  5  i 

=  f  1.  5f 

-  3,f"VMK 

=  -5i"^  K“ 

-  5  .?"^a)%K 

=  -  3.  5  f 

+  3.5i“^a;2MK 

=  +  1. 

+  1.5i”VMK 

=  + 

=  -  2.  5jC"^K^ 

=  -  1.  5i“^K^ 

=  -  2.  5i 

-2.5i 

-1  9 

=  +3. 5i 


4.  Fixed- Fixed  Beam  Sample  Solution 


Equation  (26)  represents  the  solution  of  a  uniform  symmetrical  fixed-fixed  ended 
beam.  The  square  matrix  equals  zero  for  any  natural  frequency  that  can  be  calculated 
from  the  mathematical  model.  First  the  matrix  has  to  be  expanded  to  a  single  equation. 
By  using  equation  (29)  as  a  model,  Table  I  is  constructed.  The  summation  of  Table  I 
set  to  zero  is 

OJ^  -  8u)2  M  K  +  3K2  =  0. 

If  we  use  the  quadratic  solution  equation  of 


2a 


and  set 

a  = 

b  =  -  8  M  K 
c  =  3 

then  can  be  solved  as  follows: 
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The  error  in  using  the  mathematical  model  of  equation  (26)  would  be 


Error  =  —  33'g^^^  x  100  =  2. 1%, 


A 

when  cj  is  considered.  Actually  the  error  in  terms  of  co  would  be 


^19.85-19.65 
^  19.65  19.65 


Since  the  second  value  of  as  calculated  is  about  20  times  the  first,  it  is  obvious  that 
it  does  not  represent  the  third  mode.  We  can  conclude  that  the  four  by  four  matrix  can 
be  good  only  for  the  one  lowest  natural  frequency  calculation.  If  an  eight  by  eight  ma¬ 
trix  were  used  the  accuracy  of  the  first  mode  would  be  increased  and  the  third  mode 
made  available.  Figure  13  is  the  diagram  for  equation  ( 26).  It  is  noticeable  that  only 
two  masses  are  in  oscillation.  This  is  not  sufficient  freedom  to  simulate  a  good  third 


FIGURE  13.  FOLDED  MODEL  FOR  FIXED-FIXED  BEAM 

mode.  That  is  why  the  third  mode  can  not  be  calculated  by  this  small  size  matrix. 

It  should  be  remembered  that  the  symmetrical  method  eliminates  the  second  mode. 

For  the  first  mode  the  column  matrix  of  equation  ( 26)  can  not  be  determined  in 

absolute  values.  A  solution  can  be  performed  in  the  following  manner  by  using  equation 
(26): 


20 


Substitute  in  the  calculated  value  of 

1 

Divide  the  column  matrix  by  y  . 

3 


Solve  the  set  of  simultaneous  equations  for  R(/y3,  M-(/y3  and  y/ys- 

Note  that  E/ys  is  the  reaction  force  for  the  first  mode  at  the  beam  attachment  location 
per  unit  deflection  of  the  beam  at  ys  and  that^c/ys  is  the  moment  resulting  from  the 
first  mode  at  the  beam  attachment  per  unit  deflection  at 

Other  aspects  of  the  system  can  be  calculated.  The  average  slope  of  the  beam 
for  the  first  mode  between  y2  and  y^  is 


or 


in  radians  per  unit  deflection  at  ys-  The  maximum  stress  in  the  beam  is 


S  =  -^  ,  (30) 

where  C  is  the  distance  from  the  neutral  axis  to  the  extreme  fiber.  Since  we  have  cal¬ 
culated  ^/ys  we  get 


_S 

ya 


Ys  I 


(31) 


as  the  maximum  stress  per  unit  deflection  of  y3.  It  must  be  remembered  that  stress 
calculated  by  equation  (31)  is  for  one  mode  of  vibration.  The  total  normalized  stress 
resulting  from  all  modes  would  be  the  complex  summation  of 


where  S^/y  is  the  normalized  maximum  stress  of  the  ntn  mode. 

Hinged- Hinged  Beam  Sample  Solution 

Equation  27  is  a  matrix  solution  for  a  uniform  symmetrical  beam  with  both  ends 
hinged,  or  simply  supported.  Equation  27  is  solved  using  equation  29  as  a  model.  Table 
n  is  a  tabulation  of  the  calculations. 
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TABLE  II.  EQUATION  (21)  MATRIX  EXPANSION 


AFKQ  =  .  5(K  -  w^M)  (K  -  co^M)  (-K  -  .  5a;®M^  +.  +.  Sw 


AGLN  =  .  5(-2K)  (-K)  (-Soj^M)  -3.0 

AHJP  =  .  5(K)  (-2u)%)  (K -20)%)  +2.0  -1.0 

-AFPL= -.  5(K  -  0)%)  (K  -  2a;%)  (-K)  +1.0  -1.5 

-AJGQ  = -.  5(-2o)%)(-2K)(-K  -  0)%)  +2.0  +2.0 

-ANKH  = -.5(-3o)%)(K  -  o)%)  (K)  -1.5  +1.5 

-BEKQ  =  3K(1.5)(K  -  a;%)(-K  -  o)%)  4.5 

-BGLM  =  3K(-2K)  (-K)  (3.  5) 


-BHIP=  3K(K)(2.5((K  -  2o)%)  -15.0 

+BEPL= -3K(1.  5)(K -2a)%)(-K)  -9.0 

+BIGQ  =  -3K(  2.  5)  ( -2K)  ( -K  -  oj%)  - 15. 

+BMKH  = -3K(3.  5)  (K  -  a)%)(K)  +10.5 

+CEJQ  =  K(  1. 5)  (-2o)^M)  (-K  -  o)^M)  +3.0  +3.0 

+CFLM  =  K(K  -  co%)  ( -K)  ( 3.  5)  +3.  5 

+CHIN  =  K(K)(2.5)  (-3o)%)  -7.5 

-CENL  =  -K(1.5)(-3o)%)(-K)  -4.5 

-CIFQ  =  -K(2.  5)  (K  -  a)%)  ( -K  -  oj%)  -2.  5 

-CMJH  =  -K(3.  5)(-2oj%)  (K)  +7.0 

The  summation  of  Table  II  yields; 


0  =  .  5o)%^  +  9o)^M^K  -  28.  5o)^MK^  +  2K^ 


MK^  -.  5K^ 

+.  5 

-4.  5 

+21.0 
+7.  5 
+4.  5 
-15.0 
-10.5 

-3.  5 

+2.  5 
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If  we  set 


X  K 
M 


the  above  equation  becomes 


0  =  .  5  X®  +  9  -  28.  5  X  +  2. 

We  know  that  x  will  be  a  small  number  in  the  order  of  0.  1.  Therefore  the  x*  term  in  the 
equation  can  be  neglected.  Use  the  quadratic  equation  general  solution  to  get 

_  28.5  (28.5)^  -  72 

then 


18 


0718 


2  .0718K  .0718  El 


M 


M 


The  equation  for  a  uniform  symmetrical  hinged  beam  is 


9  ■ 0752  E  1 

M 


The  error  for  is 


100  X 


0752  -  . 0718 
.  0752 


4.  5%. 


The  error  for  w  is 


100  X 


-s/Tsi' - 'sf718 

'vTtsI 


2.  3%. 


As  stated  in  the  fixed^fixed  beam  example,  the  4  by  4  matrix  is  in  reality  too 
small  for  good  problem  solution.  The  third  mode  solution  is  not  available.  The  normal¬ 
ized  moment,  end  resultant  force  and  deflections  for  the  first  mode  can  be  calculated  as 
well  as  the  average  slope  between  hinges. 

TABLE  III.  MATRIX  EXPANSION  OF  EQUATION  ( 28) 

2  K(K-co2  M)  (K-w^  M)  (-K-co^M)  =  -2K^  +  2K®to2  M  +  2K^  -  2Ka)®M^ 

2  K(-2K)  (-K)  (-Scli^M)  =  -12K®aj2M 

2K(-2a)^M)  (K -20)2  M)  =  -  4K*  0)2  M  +  8 
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Table  III  (Cont'd) 


-2  K(-3aj2  M)  (K  -  M)  (K) 

-2  K(K-w2  M)  (K  -2a;2M)  (-K) 

-2  K(-2  0)2  M)  (-2  K)  (-K  M) 


6  k2  0)2  M  -  6  k2  0)^ 

2  -  6  k2  0)2  M  +  4  k2  w*  m2 

+  8  K®  0)2  M  +  8  k2  0)^  m2 

-  6  K®  0)2  M  +  162  k2  o)^  m2  -  2  K  O)®  M® 


Free-Free  Beam  Sample  Solution 


Equation  ( 28)  is  the  matrix  solution  for  a  uniform  symmetrical  beam  with  both 
ends  free.  The  square  matrix  of  equation  (28)  reduces  to: 


2  K 


K  -  0)2  M 

-  2o)2  M 

-  3  0)2  M 


-  2  K 

+  K  -  0)2  M 
+  K  -2o)2  M 


+  K 

-  K 

-  K  -  0)2  M 


The  solution  of  this  matrix  when  set  equal  to  zero  is  from  Table  HI  and  is 
0  =  -  2  o)®  M®  K  +  16  o)^  m2  k2  -  6  0)2  M  K®. 

X  K 

If  we  let  0)2  =  — ^  ,  then  the  above  equation  reduces  to 
0  =  +  2  x2  -  16  x+6=.  395. 

Then, 

9  _  ■  395  K  _  .  395  E  I 

M  i®M  ■ 

The  first  natural  frequency  of  the  free-free  beam  should  be 


9  .386  El 

—  - - - 

1  Mi® 


The  error  in  terms  of  o)2  will  be 
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100  X 


.  395  -  ■ 386 
.  395 


2.  3%. 


The  error  in  terms  of  cu  will  be 


100  X 


\/395V386 

'/396 


1. 1%. 


Cantilever  Beam  Sample  Solution 


The  cantilever  beam  is  a  special  case  of  the  beam  equation.  It  has  one  end  fixed 
and  one  end  free.  Thus  it  is  not  a  symmetrical  beam.  Figure  14  represents  a  diagra- 
matic  model  that  fits  a  4  by  4  matrix.  This  particular  model  is  used  because  it  utilizes 
the  same  model  configurations  of  equations  26  through  28.  Moments  are  taken  at  the 
hinge  locations.  Equation  ( 32) represents  the  summation  of  moments  about  the  first 
mass,  Ml,  by  using  the  left  input  to  balance  the  moment  of  spring  K^.  Equation  (33)  lists 
the  summation  of  moments  that  are  to  the  left  of  mass  M2  and  reacting  on  spring  K2. 
Equation  (34)  represents  the  summation  of  the  moments  acting  to  the  right  of  mass  M2 
and  acting  on  spring  K2.  This  switching  of  sides  is  convenient  because  spring  K^,  ac¬ 
cording  to  this  model,  has  no  deflection.  Equation  (35)  is  the  summation  of  forces  acting 
in  the  vertical  direction.  Equations  36  and  37  equate  the  spring  deflections  to  hinge  dis¬ 
placements.  Equation  38  follows  from  equations  (32)  through  (38).  Table  IV  is  the 
solution  of  the  square  matrix  of  equation  (38). 


A  A2  A3 


FIGURE  14.  MODEL  OF  A  CANTILEVER  BEAM 
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Cantilever  equations 


-A1K1+— -  .5Ro 


-A2K2  +  -P-  1.5  Ro 


-A2K2  +  w  M3y3 

-Ro  +  W%2y2  + 


Ai  ="  y2 

A2  =  YS  -  2y. 


•K  +  cj2m  y2 


TABLE  IV.  MATRIX  EXPANSION  OF  EQUATION  (38) 


(-1.5)  (2K)  (OJ^M) 


-3  OJ^K 


(2K)  (-1)  (-K  +  oj^M) 


2K^i 


-2  aj2  K  M  i 


ri  (0) 


■i-V(-l)  (2K)  (-K) 


-2  i 


(-1.5)  (w^M)  (-K  +  co^M) 


-1.5a;2KMi“^  +1.5aj^M^i"^ 


-i-^  (0) 


.5  (i~^)  (+2  K)  (0)2  M) 


+  0)2  K  M  i 


.5  (0) 
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Table  IV  (Cont'd) 


.5  (0) 

-.5  (0) 

-.5  (0) 

5  M)  (-K  +  (J-  M) 

-K  (0) 

-K  (0) 

-K  (0) 

+K  (0) 

K  (0) 

K  (-1)  (-K  +  M) 

Total  of  Table  IV 
Setting  the  square  matrix  equal  to 
0  =  -  6  cc)2  K  M  +  kI 


0 

0 

0 

+0.5a;2KMi"^  -0. 

0 
0 
0 
0 
0 

-0)2  KMi”^ 

K^i"^  -ecr^KMi"^  +w 
we  get 


Let 


CO 


2  - 


X  K 
M 


then 


0  =  x^”6x+l, 


X  = 


6  ±  \/  36  -  4 
2 


6  ±4^ 
2 


6  -  5.  67  .33 

2  “2 


.  165, 


then 


m2 
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,  .  165  K  .  165  E  I 
1  M  M 


The  standard  solution  for  the  first  mode  of  a  cantilever  beam  will  give 
,  . 153  E  I 

=  - 5— 

1  M 

The  error  in  terms  of  is 


.  165  -  .  153 
.  153 


X  100  =  7.  9%  . 


The  error  in  terms  of  clj  is 


163  -  '\/l53 
's/153 


X  100  =  3.2%. 


Since  there  are  only  two  displacement  values  in  this  example,  y2  and  y3,  good  accuracy 
can  not  be  expected. 


Summary  of  the  Worked  Examples 

The  four  examples  worked  here  are  algebraic  with  the  end  resultant  in  a  form 
that  could  be  compared  to  the  general  beam  natural  frequency  equation.  For  an  example 
where  the  masses,  stiffnesses,  and  length  dimensions  can  be  determined,  the  Eigenvalue 
solution  of  the  square  matrix  will  find  the  natural  frequencies  directly  without  the  substi- 
E  I 

tution  of  “p-  in  for  K  at  the  end  of  the  problem.  Actually  this  substitution  will  be  made 

when  determining  the  initial  values  of  stiffness  to  insert  into  the  matrix.  It  need  only  be 
made  once. 

Although  a  4  by  4  matrix  is  obviously  too  small  to  assure  enough  input  detail, 
the  accuracy  of  the  four  problems  is  good  enough  to  point  out  the  validity  of  the  method. 


Non-Uniform  Beam  Solutions 

The  mathematical  model  put  forth  in  this  paper  is  adaptable  to  beams  of  non- 
uniform  stiffness  and/ or  mass  as  well  as  beams  with  a  curved  neutral  axis.  These  prob¬ 
lems  can  be  developed  to  fit  a  staudard  Eigenvalue  solution  to  a  square  matrix  but  the 
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matrix  itself  is,  of  course,  not  standard  in  form.  Each  of  these  special  situations  must 
have  a  special  set  of  equations  developed.  First  we  must  examine  the  limitations  of  the 
method.  For  an  example,  consider  Figure  15  as  typical  of  a  special  beam  base. 


FIGURE  15.  NON-UNIFORM,  NON  SYMMETRICAL  BEAM 

Each  link  of  the  simulated  model,  Figure  16,  must  have  like  dimensions  but  the  indivi¬ 
dual  values  of  M  and  K  can  vary.  The  trussed  links  must  be  small  enough  so  that  the 
net  deflection  of  the  beam  hinge  can  be  considered  a  net  deflection  of  the  spring.  Error 
will  be  proportionate  to  the  relative  angle  mode  by  the  bases  of  adjacent  trusses.  This 
angle  will  increase  as  the  length  "1"  is  increased,  and  decrease  as  the  length  "i"  is  de¬ 
creased.  The  allowable  error  puts  a  limitation  on  the  model  truss  size  and  therefore 
on  the  minimum  number  of  trusses. 


FIGURE  16.  NON-UNIFORM  BEAM  MODEL 

The  deflection  at  Mj  will  have  a  resultant  with  both  an  "x"  and  a  "y”  component. 
However,  if  a  single  coordinate  system  of  response  relative  to  location  is  worked  out, 
a  smaller  matrix  results. 


29 


In  this  particular  model  there  is  no  moment  on  spring  K^.  The  motion  of  each 
mass  concentrated  at  the  hinge  is  considered  to  be  in  the  direction  of  a  ray  that  bisects 
the  angle  between  the  adjacent  truss  bases.  The  moment  arm  about  will  then  be  a 
perpendicular  struck  through  the  location  of  to  the  directional  ray.  By  using  this 
method  a  set  of  moment  equations  can  be  developed  in  the  same  manner  as  for  the  straight 
beam.  The  force  summation  must  equal  zero  in  both  the  x  and  the  y  direction.  A  set  of 
equations  can  be  developed  to  set  up  a  matrix  equation  that  will  solve  for  Eigenvalues  and 
Eigenvectors. 

The  use  of  the  mathematical  model  presented  herein  probably  has  its  greatest 
value  in  helping  to  solve  the  more  complicated  situations  by  reducing  the  specifics  to 
a  model  that  can  be  solved. 


SECTION  IV.  VARIATIONS  ON  THE  ANALYTICAL  BEAM  MODEL 


Initial  Displacement 

If  an  initial  displacement  is  known,  it  replaces  the  appropriate  y.  For  example: 
if  y4  has  an  initial  displacement  of  .  005  inches,  this  value  will  replace  y4  in  the  column 
matrix  of  equation  (23) . 


Mode  Shape 

The  maximum  displacement  is  given  an  initial  displacement  of  unity.  A  plot 
of  the  values  of  y  versus  location  on  the  beam  will  describe  the  mode  shape.  Curve 
fitting  procedure  will  equate  y  to  x. 

Forced  Oscillation 

Forced  oscillations  will  enter  the  matrix  equation  (23)  as  a  column  matrix  on  the 
left  side.  The  procedure  is  illustrated  by  Figure  17. 


FIGURE  17.  FORCED  EXCITATION  EXAMPLE 
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0 


+  .  5Fi 
+  1.5  Fi 
+  2.  5  Fi 
+  .5  F2 
+  (Fj  +  F2) 


Mr 

=  -  AiKi  +  -f  -  .  5Ro 

=  -  A2K2+—  -  1.5 
M- 

=  - A3K3+-p-2.5R„  +  cd2  M2y2 

M- 

=  -  A4K4  +  ^  -  3.  5  Rfl  +  20)2  +  ^2 

=  *“  A4K4 

=  -  Rq  +  0)2  M2y2  +  1^2  M^2  ^sYs  +  M4y^ 


Ai-y2 
A2  =  Yg  -  2Y2 
As  =  Y4  -  2Ys  +  Y2 
A4  =2y^-  3Y4  +  Y3 


0 

-  .5 

+i-l 

-Ki 

0 

0 

0 

R 

0 

+  .  5F 4 

-1.5 

+i-l 

+2K2 

-K2 

0 

0 

+  1,  5Fi 

-2.  5 

+i-l 

(-K3  +  0)2  M2) 

+2K3 

-K3 

0 

^2 

+  2.  5F4 

-3.  5 

+r^ 

+  2o)2  M 

(-K4  +  o)2m) 

+  3K4 

-2K4 

y3 

+  .  5F2 

0 

0 

0 

-K4 

+3K4 

-2K4 

y4 

+  (Fi  +  F2) 

-  1 

0 

0)2M2 

0)^M3 

0)2M4 

0 

^5 

Mixed  End  Conditions 

When  the  end  of  the  beam  is  anchored  to  a  flexible  member,  i  Aq  Kq  is  sub¬ 
stituted  in  equation  (23)  for  Mq  and  Ao  =  2y^;  is  then  a  variable  and  replaces ^0  “ 
right  side  column  matrix. 
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When  the  end  condition  has  a  rotational  oscillation,  f(XQ)  is  then  a  known  input 
f(x^)  =  2yi  and  its  function  is  substituted  into  the  right  hand  column  matrix  for  y^. 

When  the  end  is  oscillating  in  the  y  direction,  f(yo)  is  a  known  function  and  is 
substituted  into  the  right  hand  column  matrix  (Fig.  18) . 


Multi-Supported  Beams 


For  a  multi  supported  beam,  Rj^  must  eliminate  in  the  matrix  equation.  Thus 
Rjj  will  replace  co^  yj^  in  the  basic  equations  and  R^^  will  replace  y^  in  the  right  hand 
column  matrix  as  an  unknown  (Fig.  19) . 


FIGURE  19.  EXAMPLE  OF  MULTI  SUPPORTED  BEAM 
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In  the  non-rigid  supported  beam,  as  illustrated,  the  support  must  be  at  a  hinge 
between  triangles.  In  each  instance  of  support,  for  example  where  ct>^  Mj  and 
cjo^  M5  yj  are  in  the  matrix  (23) ,  the  quantities  (co^  -  kJ)  and  M5  -  Kg)  are  re¬ 
spectively  substituted  since  the  spring  will  resist  motion  at  the  support.  If  the  springs 
are  initially  not  at  rest,  the  respective  forces  must  be  inserted  (Fig.  20) . 


K  K  K  K  K  K  K 


FIGURE  20.  EXAMPLE  OF  FLEXIBLE  SUPPORTED  BEAM 

SECTION  V.  ADAPTION  OF  A  BEAM  ANALYTICAL  MODEL 
TO  A  PLATE  ANALYTICAL  MODEL 


The  beam  of  Figure  3  is  simulated  by  a  model  that  is  most  convenient  for  all  of 
the  beam  problems.  Other  models  may  be  constructed  to  meet  specific  purposes.  The 
beam  of  Figure  3  is  not  readily  adaptable  to  plate  configurations.  Therefore,  it  is 
necessary  to  modify  our  approach  and  construct  a  basically  new  model  as  shown  in 
Figure  21.  This  new  model  represents  a  change  in  mass  attachment  that  complicates 
the  coupling  of  mass  to  the  column  matrix  of  variables. 

Kj  K2  Kg  K4 
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The  displacements  of  Mi,M2  and  M3  are  the  average  of  the  two  adjacent  hinge 
displacements.  Thus  the  dynamic  force  of  may  be  represented  by 

In  order  to  expand  the  beam  to  a  two  dimensional  model,  the  triangle  concept  is 
replaced  with  a  massless  rectangular  pyramid  as  shown  in  Figure  22  and  Figure  23. 


FIGURE  22.  PYRAMID  CONFIGURATION 

In  the  development  of  the  matrix  equation  for  this  beam  the  mass  displacement 
may  be  taken  as  average  of  two  diagonal  adjacent  hinges;  for  convenience,  the  top  left 
and  bottom  right  of  each  pyramid  base  are  taken.  The  displacement  of  spring  is 


and  the  displacement  of  is 


The  first  six  lines  of  matrix  equation  (39)  is  the  analytical  model  of  the  two- 
dimensional,  four-spring  beam.  This  model  has  no  advantages  for  beam  problems  over 
equation  23  and  is  more  complicated;  but  the  corner  hinges  allow  for  modal  coupling  of 
a  plate  problem. 
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Equation  (39) 


>>>5  >l>i  >>  >>  >>>^ 
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These  pyramids  are  linked  at  the  corners  to  form  a  beam  as  in  Figure  23 


TOP  VIEW 


K2  K3  K4 


SIDE  VIEW  FOR 
NOMENCLATURE 


TOP  VIEW  FOR 
iNOMENCLATURE 


FIGURE  23.  PYRAMIDIC  MODEL  SUITABLE  FOR  PLATE  ADAPTION 


The  plate  made  up  of  four- spring  beams  may  be  represented  as  in  Figure  24. 


FIGURE  24.  FOUR  SPRING  PLATE  MODEL 

There  are  three  horizontal  beams  and  three  vertical  beams.  The  four  corners  of  the 
plate  are  not  a  part  of  any  beam.  However,  when  the  assumption  that  the  pyramids  are 
rigid  is  applied,  the  three  hinges  of  the  corner  segments  that  are  attached  to  beams  will 
determine  the  displacement  of  the  corners. 

Because  of  the  coupling,  the  number  of  displacement  variables  of  the  plate  is 
not  equal  to  the  sum  of  the  six  beams  and  the  four  corners.  This  fact  reduces  the 
number  of  required  equations  to  make  up  a  plate  matrix  equation.  Each  beam  model 
must  include  all  of  the  variables  of  that  beam.  When  fewer  equations  are  required,  the 
fifth  and  sixth  lines  of  the  equation  may  be  deleted  and  line  seven  substituted  as  the  fifth 
equation  of  the  beam  set.  The  plate  of  Figure  24  will  require  a  rectangular  matrix  of 
60  by  36  which  reduces  to  36  by  36  when  the  end  conditions  are  inserted.  Equation  40 
gives  the  analytical  model  of  the  plate.  The  end  conditions  are  selected  in  the  same 
manner  as  the  single  beams  and  all  the  variations  as  described  for  the  beams  may  be 
applied  to  the  plate  matrix. 

The  stiffness  constant  of  the  two  dimensional  beam  model  is 
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where  p  is  Poisscm's  Ratio  and  the  other  nomenclature  is  the  same  as  in  the  beam. 

The  4  by  4  spring  plate  solutions  will  be  accurate  for  only  the  first  mode  of  the 
plate.  Practical  considerations  would  dictate  a  matrix  up  to  the  capacity  of  the  available 
computer.  For  example  a  plate  with  13  by  13  springs  and  good  for  3  modes  in  each 
direction  would  require  a  169  by  169  matrix;  a  plate  good  for  4  modes  in  each  direction 
would  require  a  289  by  289  matrix. 

Since  maximum  amplitudes  will  be  encountered  in  the  lower  modes,  there  is 
practical  value  from  the  use  of  the  model.  The  Eigenvalue  problems  need  only  to  be 
computed  for  a  limited  number  of  the  lowest  frequencies.  Higher  frequency  solutions 
would  have  little  meaning. 


SECTION  VI.  VARIATIONS  OF  THE  PLATE  MODEL 


There  are  two  major  variations  of  the  plate  model  that  have  significance.  There 
are  adaptations  of  the  model  to  allow  the  inclusion  of  stiffening  beams,  and  the  simula¬ 
tion  of  a  cut-out  plate  for  a  complex  beam  structure.  Both  represent  a  use  of  additions 
or  subtractions  of  springs  and  masses  to  the  plate  model.  The  matrix  equation  can  be 
modified  to  meet  the  specific  applications. 

Shell  Model 


The  analytical  model  can  be  modified  to  simulate  a  thin  cylindrical  shell  struc¬ 
ture.  The  deflection  is  assumed  in  a  radial  direction.  With  r  0  substituted  into  the 
matrix  equation  for  i  of  the  wraparound  beam,  the  wraparound  beams  are  continuous 
and  therefore  have  no  end  conditions.  This  changes  the  number  of  unknowns  and  also  the 
character  of  available  equations  to  set  up  a  suitable  matrix. 

The  adaptation  of  analytical  principles  used  in  this  paper  will  be  extended  to 
cylinders  and  published  at  a  later  date. 
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SECTION  VII.  CONCLUSIONS 


The  dynamics  of  a  beam  or  plate  can  be  simulated  by  a  model  made  up  of  mass¬ 
less  constraints,  springs,  and  concentrated  masses.  By  unique  arrangement  of  the 
mechanical  elements,  a  mathematical  model  may  be  produced  for  structures  not  gener¬ 
ally  solvable  by  usual  techniques.  The  analytical  approach  has  been  derived  and  the  re¬ 
sults  proven  by  the  use  of  simple  hand-solvable  models.  Extension  to  more  complicated 
structures  is  possible  using  digital  computers  for  matrix  Eigenvalue  and  Eigenvector 
solutions. 

George  C.  Marshall  Space  Flight  Center 

National  Aeronatuics  and  Space  Administration 
Huntsville,  Alabama,  February  1965 
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